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Abstract 

We study the regularities in the multiparticle production data obtained from different 
types of collisions indicating the universality of the hadroproduction process. The sim- 
ilarities of such bulk variables like the charged particle mean multiplicity and the pseu- 
dorapidity density at midrapidity measured in nucleus- nucleus, (anti)proton-proton 
and e + e~ interactions are analysed according to the dissipating energy of participants 
and their types. This approach shows a good agreement with the measurements in a 
wide range of nuclear collision energies from AGS to RHIC. The predictions up to the 
LHC energies are made and compared to experimental extrapolations. 



1. Nucleus-nucleus collisions at RHIC probe matter at high densities and temperatures 
ever reached and provide us with an opportunity to investigate strong interactions up to ex- 
tremely high energy density parton collisions. Of fundamental interest are bulk observables 
such as multiplicity and particle densities (spectra), which are sensitive tools for probing 
the dynamics of strong interactions. They give us information about the system formed in 
high energy collisions, both after cooling and hadronisation and as well during the formation 
and evolution of the collision initial state, and are thus powerful in distinguishing between 
different particle production models. Recent measurements at RHIC revealed striking evi- 
dences in the hadron production process including some universality between such basic ob- 
servables as the mean multiplicity and the midrapidity density in complex ultra-relativistic 
nucleus-nucleus collisions vs. those measured in relatively "elementary" e + e~ interactions. 
The values of these bulk observables are found to be similar for both types of reactions when 
measurements are normalised to pairs of participants ( "wounded" nucleons PQ in heavy ion 
collisions) at the nucleon-nucleon center-of-mass (cm.) energy, ^/snn, comparable to the 
cm. energy ^/s^T of e + e~ annihilation 012] • This phenomenon is found to be independent 
of the energy spanning from ^/snn = 19.6 GeV to 200 GeV at RHIC. Assuming a universal 
mechanism of hadron production in both types of interactions which then depends only on 
the amount of energy transformed into particles produced, one would expect the same value 
of the observables to be obtained in hadron-hadron collisions at close cm. energies. How- 
ever, this is not the case: comparing measurements in hadronic data jHEj to the findings 
at RHIC, one obtains 00 0111 El 0110111111121 113 El quite lower values in hadron-hadron 
collisions. In the meantime, recent RHIC data from deuteron-gold interactions at ^snn 
21)0 GeV unambiguously points to the same values of the mean multiplicity as measured in 
antiproton-proton collisions 00111111^11101 • 

In this paper, we give a phenomenological interpretation of the above regularities obtained 
and compare these to the earlier data at lower energies. The particle production is considered 
to depend on the amount of the participant energy dissipated in collision. The comparison 
between different reactions is made on the basis of the type of participating patterns. This 
consideration introduces a tripling factor to be taken into account for the cm. energy of a 
reaction and for appropriate treatment of participants to correctly compare the bulk variables 
from different reactions. Good agreement between the description proposed here and the 
experimental measurements is found. 

2. In the consideration given here, the whole process of a collision is interpreted as 
the expansion and break-up into particles of an initial state, in which the whole available 
energy is assumed to be concentrated in a small Lorentz-contracted volume. There are no 
any restrictions due to the conservation of quantum numbers besides energy and momentum 
constraints allowing therefore to link the amount of energy deposited in the collision zone 
and features of bulk variables in different reactions. This approach resembles the Landau 
phenomenological hydrodynamical description of multiparticle production in relativistic par- 
ticle collisions ^7j. Though the hydrodynamical description does not match ideally the data 
on multiparticle production in the whole range of pseudorapidity, the model-based energy- 
scaling law of multiplicity gives good agreement with the measurements in such different 
reactions as nucleus-nucleus, pp, e + e~ and up collisions [T%1IT9*] . Recently, the Landau model 
prediction for the longitudinal pion transport has been observed to be well within 5% ac- 
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curacy with the heavy-ion data [20] • This indicates that the main asserts of the Landau 
approach are useful to estimate fractions of the energy dissipated into particles produced in 
different reactions, particularly in nucleus-nucleus collisions |2*T] . 

As soon as the collision of two Lorentz-contracted particles leads to the full thermali- 
sation of the system before extension, one can assume that the production of secondaries 
is defined by the fraction of energy of participants deposited in the volume of thermalized 
system at the collision moment. This implies that there is a difference between results of 
collisions of structureless particles like electron and positron and composite particles like 
proton, the latter considered to be built of constituents. Indeed, in composite particle colli- 
sions not all the constituents deposit their energy when they form a small Lorentz-contracted 
volume of thermalised initial state. As a result, the leading particles [22] formed out of those 
constituents which are not trapped in the interaction volume, carry away a part of energy 
effectively making it unavailable to participate in production of secondaries. From the other 
side, colliding structureless particles are ultimately stopped as a whole in the initial state 
of thermalized collision zone depositing their total energy in the small Lorentz-contracted 
volume. The latter makes this energy of the incoming particles wholly available for the 
production of secondaries. 

A single nucleon represents a superposition of three constituent quarks. In this additive 
quark picture [2H1, most often only one quark from each nucleon contributes to the interac- 
tion, while other quarks are considered to be spectators. Thus the initial thermalized state 
which is responsible for the number of produced secondary particles is pumped in only by 
the energy of the interacting single quark pair. The quark spectators form hadrons being not 
stopped in the thermalized volume at the collision moment do not participate in secondary 
particle production. The later means that only 1/3 of the entire nucleon energy is available 
for particle production. In e + e~ annihilation, the incident particles are structureless and 
therefore, the total interaction energy is deposited in the thermalized collision zone. From 
this and the above consideration, one expects that the resulting bulk variables like the mul- 
tiplicity and rapidity distributions should show identical features in proton-proton collisions 
at the cm. energy -/s pp and e + e~ interactions at the cm. energy -y/s^T ~ ^/s pp /3. Note that 
for the mean multiplicity, similar behaviour was obtained experimentally in the beginning 
of LEP activity [21]. 

In heavy ion collisions, more than one quark per nucleon interacts due to the large size of 
nucleus and the long path of interactions inside the nucleus. The more central the nucleus- 
nucleus collision is, the more interactions occur and the larger energy amount is spent for 
secondary particle production. In central nuclear collisions, a contribution of constituent 
quarks rather than participating nucleons seem to determine the properties of produced 
particle distributions [23]. In the most central collisions, the density of matter is so high 
(almost saturated) that all three constituent quarks from each nucleon may participate nearly 
simultaneously in collision depositing their energy coherently into thermalized collision zone. 
The total entire energy of nucleons included in the most central fraction of colliding nuclei 
is available for secondary particle production. Recalling that in proton-proton collisions, 
where only one out of three constituent quarks from each proton interacts, one expects the 
features of the bulk variables per pair of participants measured in the most central heavy-ion 
interactions at the cm. energy t/s^n to be similar to those from proton-proton collisions at 
three times larger cm. energy, A /s^ p " ~ 3 v /5nn. This makes the most central collisions of 
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nuclei akin to e + e collisions at the same cm. energy from the point of view of the resulting 
bulk variables. 

3. Applying the above consideration, we compare in Fig.Q in the a/snn energy range from 
a few GeV to 200 GeV, the cm. energy dependence of the mean multiplicity as measured 
in nucleus-nucleus 00121112111211121, and e+e" interactions Hp EH EH1 [13 101 EM 131111 
vs. that but at the three times larger cm. energy in (anti)proton-proton pp/pp collisions 
|30p3ip32p33p34j . The multiplicity values in nuclear collisions are given divided by the number 
of participant pairs and the energy value gives the cm. energy per nucleon, a /snn. For 
the data on e + e~ annihilation at energies above the Z° peak, we give the multiplicity values 
averaged 1 here from those recently published by LEP experiments at LEP1.5 y/s^ = 130 GeV 
jSllEIl and LEP2 v /s^ = 200 GeV (HUGH energies: 23.35 ± 0.20 ± 0.10 (LEP1.5) and 
27.62 ± 0.11 ± 0.16 (LEP2). Figure shows also the average multiplicity fit to pp/pp data 
obtained in an d the ALEPH fit [37], based on the 3NLO perturbative QCD calculations 
[IB], to e + e~ data. 

First, one can see from Fig. Uthat the pp/pp data are very close to the data from e + e~ 
annihilation measured at the cm energy t/s^ = y/s pp /3. This nearness, as the fits confirm 
spanning from GeV to TeV cm. energies, decreases the already small deficit in the e + e~ 
data compared to the pp/pp data as the energy increases. Note that at TeV energies this 
deviation is almost five-six times smaller the area covered by the fit errors. The deviation 
can be attributed to the inelasticity factor, or leading particle effect |22j in pp/pp collisions, 
which is known to decrease with the cm. energy. Then, at lower cm. energies, some fraction 
of the energy of spectators contributes more into the formation of the initial state, while the 
spectators passed by. This leads to the excess of the average multiplicity in pp/pp collisions 
compared to that in e + e~ interactions at ^/s^ = y/s pp /3 as it is observed from Fig. [TJ 

Comparing further the average multiplicity values from pp/pp collisions to those from 
heavy-ion reactions, one finds that the data points are amazingly close to each other when 
the heavy-ion data is confronted with the hadronic data at the energy y/s pp = 3 ^/snn and, 
according to the just made observation, to the e + e~ data at the same energy as the nucleus- 
nucleus reaction data are taken. All these findings agree well with our interpretation of the 
multihadron production process and the consequences the bulk variables to be similar in the 
reactions considered when measured at ^/s~^ and ^ s /s^ being about a third of y/s pp . The 
inclusion of the tripling energy factor indeed allows energy-independent description of such 
a fundamental variable as the mean multiplicity in simultaneously e + e~, pp/pp and central 
nucleus-nucleus collisions. This shows that the particle production process in headon nucleus- 
nucleus collisions is derived by the energy deposited in the Lorentz-contracted volume by 
a single pair of effectively structureless nucleons similar to that in e + e~ annihilation and 
of quark pair interaction in pp/pp collisions. To note is that an examination of Fig. ^ 
reveals that there is no need to rescale the y/s^ by a factor of 1/2 to match the pp data 
as earlier was assumed [2111111 , while recognised [Hj to unreasonably shift the e + e~ data 2 on 
the pseudorapidity density at midrapidity when compared to the heavy-ion measurements. 

1 The averaging procedure is adopted from the PDG review 46] for the averaging over N measurements. 
The errors are multiplyed by the S-f actor, where S = \/x 2 /(N — 1) > 1 or S = 1 otherwise. 

2 Recall that the factor 1/3 but not 1/2 was already found earlier in [21] to give an appropriate rescaling 
of the pp mean multiplicity data relative to those from e + e _ annihilation. 
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This discrepancy, as shown, finds its explanation in our consideration, within which the data 
on both the mean multiplicity and the midrapidity density (vide infra) are self-consistently 
matched for different reactions. 

Fig. ^ shows that the mean multiplicities measured in different types of interactions are 
close to each other starting from the SPS a/snn energies, and become particularly close 
at v /Inn ^ 50 GeV. At lower energies, however, small deviation in the measurements is 
visible. The nucleus-nucleus data points are slightly below the measurements from e + e~ and 
hadronic experiments. At snn ^ 10 GeV, the nuclear data departure increases further, 
and the data start to decrease faster with the energy decreasing than the measurements from 
pp and e + e~ experiments do. On the other hand, as cm. energy increases above a few tens 
GeV, the heavy-ion data start to exceed the e + e~ data and reaches the mean multiplicity 
values from pp interactions. From this, one concludes on the two different regions of the rise 
of the mean multiplicity nuclear data with the energy increase: one being steeper at a/snn 
of a few tens GeV, and another one slower at lower energies down to a few GeV. One also 
can see that the multiplicities in nuclear reactions increase faster with the cm. energy than 
this can be found in other interactions, while the effect is relatively small keeping the data 
points close to each other at fixed cm. energy ^snn and ^/s^. of the value of y / s^p/3. 

The observations made can be understood in terms of the overlap zone. At relatively 
low energies, the initial thermalized state, formed in headon nuclear collisions, are most 
likely to occur in the overlap zone of colliding nuclei along the incident direction, so the 
nucleons at the periphery are not contributing. Since the energy of the collision is small and, 
therefore, is almost entirely and fast converted into particles produced, the rest parts of the 
nuclei start to fragment into pieces by expanding overlap zone. The smaller the energy is, 
the less particles will be produced in a collision. At low energies, due to our approach, the 
mean multiplicity in pp collisions at ^s pp is expected to be larger than in nucleus-nucleus 
collision at a/snn = ^ s a/ s nn increases, the mean multiplicity in nucleus- nucleus 

collisions is supposed to increase faster than that in pp/pp interactions at -/s pp = 
Indeed, in pp/pp collision we consider a single pair of constituent quarks to form participants, 
while in nuclear collisions, the effectively structureless protons entirely participate. So, at 
the same energy per participant, at low energies, in (anti)proton collisions, all the energy of 
a pair of interacting constituent quarks is converted into particles produced, while in nuclear 
collision some fraction of energy is used to break-up nuclei by expansion of the small overlap 
zone. As the energy increases, the nature of the rise in (anti)proton collisions should not 
change significantly since (anti)protons are considered interacting the same way as at lower 
energies, i.e. via a pair of constituent quarks depositing their energy into secondaries. In 
contrary, in (central) nuclear collisions, the Lorentz contraction starts to play a significant 
role involving more nucleons into overlap zone. Due to this, more energy becomes available 
for particle production which makes the mean multiplicity rise faster than at lower energies. 
To add is that this effect attracts currently special consideration in different approaches pE]. 

Note that the total multiplicity is not very sensitive to the above, and this is seen in 
Fig. ^ where, as we already mentioned, the difference between hadronic and nuclear data 
is small, once tripling energy factor is taken into account. On the other hand, due to fact 
that differences are related mostly to the interaction zone and consequently, to the central 
rapidity region, one would expect the differences to be more pronounced in the midrapidity 
comparison, as we indeed find below. 



4 



4. In Fig. we compare the values of the pseudorapidity densities per participant pair 
at midrapidity as a function of cm. energy, from most central nucleus-nucleus collisions at 
RHIC jgllimE^imEni; CERN SPS |H21IMHMH^1 ; and AGS ^ vs. those measured in pp/pp 
interactions at CERN [4ll30p57|l32j and Fermilab [51l34|l59j . The comparison is given the same 
way as the mean multiplicity is shown in Fig. Q i.e. the data from nucleus-nucleus collisions 
are plotted at the energy ^/snn = i/Spp/3. One can again see that up to the existing v /snn, 
the data from hadronic and nuclear experiments are close to each other being consistent 
with our interpretation. The measurements from the two types of collisions coincide at 
8 GeV < ^/snn < 20 GeV and are of the magnitude of the spread of heavy-ion data points 
at the highest energy of 200 GeV. For a few GeV energy shown, one can also see the visible 
difference. The data shown in Fig. El indicates that the deviation between the two types of 
collisions increases with the cm. energy due to faster increase of the midrapidity density 
values obtained in heavy-ion collisions in comparison with those measured in pp interactions. 
At lower energies too, the nuclear data, being lower than the pp data, increases faster. The 
latter means that, as we discussed above, in contrast to the mean multiplicity, which is in 
general defined by the total yield of the reaction, so being less sensitive to reaction details, 
the midrapidity density depends on some additional factor. As the midrapidity density is 
measured in the very central region, where the participants longitudinal velocities are zeroed, 
it is natural to assume that this factor is related to the size of the Lorentz-contracted volume 
of the initial thermalized system determined by participating patterns. 

To take into account the corresponding correction, let us consider our picture in the 
framework of the Landau model which reasonably well describes the bulk variables measured 
and is, by its nature, near to our interpretation as discussed above. Using this model, one 
finds for the ratio of the charged particle rapidity density p(y) = (2/N pajTt )dN ch /dy per 
participant pair at the midrapidity value y — in heavy-ion reaction, pnn, to the density 
p pp in pp/pp interaction, 

Pnn(O) _ 2iV ch ALp7 

Ppp(o) iv part iv c p h p Vl nn 

Here, A part is the number of participants in heavy-ion collision, N c ^ (N^) is the multiplicity 
in nucleus- nucleus (pp/pp) collision and L = ln^, where m is the mass of a participating 
pattern, e.g. of a proton, m p , in central heavy- ion collisions. According to our interpretation, 
we compare in the ratio (JTJ) the rapidity density Pnn(0) at ^/snn to the rapidity density p pp (0) 
at ^/Spp/3. Due to the above, we consider a constituent quark of mass |m p as a participating 
pattern in pp/pp collisions, and a proton as an effectively structureless participant in most 
central nucleus-nucleus collisions. Then, from Eq. one obtains: 

(n\ ^ 2 A ch [ 41n3 

Pnn O =p pp — -pp Jl - -. 2 

Apart iv^f y ln(4m^/s N N) 

Using the fact that the transformation factor from rapidity to pseudorapidity does not 
influence the above ratio and substituting the data values of and shown in Fig.^and 
of Ppp(0) shown in Fig. one obtains from Eq. the values of pseudorapidity density in 
nucleus- nucleus collisions. These values are displayed in Fig.Elby solid line. One can see that 
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the correction made provides a good agreement between the calculated /?nn(0) values and 
the data. This justifies the above argued midrapidity multiplicity dependence on participant 
type determining the volume of the initially thermalized system. Eq. (J2J) shows also the 
importance of the correction, in particular in case of the Landau model, for the participant 
type to be introduced to properly estimate the midrapidity density. This correction is in 
agreement with our description where the proper participants have to be considered for the 
corresponding reaction. To note is that the pseudorapidity density at midrapidity in e + e~ 
annihilation is shown [3| to coincide with that from nucleus-nucleus data in the wide energy 
range and is in agreement with our assumptions and the above expectations. 

As we discuss above for the mean multiplicity, one can see that the pseudorapidity den- 
sity at midrapidity is indeed more sensitive to the differences in the multiparticle production 
process at lower and higher energies. The increase of the density measured in nuclear colli- 
sions at ^/snn less than the highest SPS energy is steeper than it is above it. This is also 
true for pp/pp data which is fitted with different log functions: with the single log fit at 
lower ^Spp and with 2nd order log polynomial at higher energies [3]. It is worth noting 
how well our consideration of tripling energy and properly treating the type of participants 
- quarks in pp/pp and protons in nuclear collisions - takes into account the different energy 
behaviour in pp/pp collisions from which the nuclear data can be deduced. Note that the 
same two regions we find here studying the pseudorapidity density at midrapidity Pnn(O), 
recently has been reported in [TIHITTj by PHENIX, when analysing the ratio of the transverse 
energy density to the pseudorapidity density at midrapidity as function of energy. From this 
finding, one can expect to obtain similarities in the transverse energy density in pp/pp vs. 
that in headon heavy-ion collisions in the frame of our description. The transition region 
at the SPS energies has been found also by NA49 |2H], which can be treated within our 
approach too, without any additional assumptions. 

From the above, we conclude that the description proposed here allows explaining the 
similarity obtained between the bulk variables measured in central heavy ion collisions com- 
pared to other types of reactions. Namely, we consider the constituent quark structure as 
a frame governing the formation of the initial thermalized state in the collision zone. Then 
a pair of structureless participating patterns being stopped in a Lorentz-contracted volume 
is the main source of energy dissipating into particles produced. This makes only a fraction 
of the total energy of colliding composite object like (anti)proton available for formation of 
secondary particles. In contrary, in headon collision of nuclei, nucleons are considered as be- 
ing structureless, similar to e + e~ annihilation, owing to almost coherent deposition of energy 
by all participating constituent quarks. This provides the total energy per nucleon being 
available for the formation of the thermalized Lorentz-contracted volume and consequently, 
for the secondary particle production. Therefore, in headon nuclear collisions, and in e + e~ 
annihilation as well, a tripling factor is needed to be taken into account for energy and mass 
of interacting patterns to adjust the global variables to those in pp/pp collisions. It is worth 
noting that recently, the constituent quark picture has been exploited to reasonably model 
the heavy- ion pseudorapidity and transverse energy data 

5. An intriguing issue, which can be considered from the interpretation given here in 
tracing similarities of the bulk observables in nucleus-nucleus collisions and in "simpler" 
reactions, particularly in such like nucleon-nucleon ones, is to make some predictions for the 



6 



bulk variables in asymmetric nucleon-nucleus collisions. In the frame of our consideration, 
in these interactions, the mean multiplicity value per pair of participant measured at some 
y/^NN is expected to have the same value as that is in pp/pp collisions at the same y/s pp as 
v /snn. Indeed, assuming an incident proton in p-nucleus collisions interacts in a way it would 
interact in pp collision, the secondary particles in the reaction are considered to be created 
out of energy deposited by interaction of a single pair of constituent quarks, one of which is 
from the proton and another one from a nucleon of the interacting nucleus. So, in fact, only 
this pair of constituent quarks converts their energy into secondary particles similar to that 
in pp/pp interactions. This also implies that the mean multiplicity being defined mostly by 
the energy deposited by participants - by the pair of constituent quarks - is not expected to 
depend on the centrality of nucleus-induced collisions, i.e. on the number N pa _ n of participants 
(within uncertainties due to intranuclear effects, e.g. Fermi motion). These expectations are 
recently shown to be well confirmed in the RHIC data on deuteron-gold interactions at 
a/snn = 200 GeV - Moreover, the effect of the similarity obtained at RHIC is 

shown j3[T31[Tn] to be true for hadron-nucleus collisions at ^/s^ ~ 10-20 GeV too. The 
same seems to be correct also for the pseudorapidity density at midrapidity, which is already 
supported to be a trend piTp] . This can also explain the necessity of the multiplication [UU] of 
the pp/pp pseudorapidity distribution by a number of the Au participants to correct the QCD 
saturation model predictions [SI] to fit the data gold nucleus region of the pseudorapidity 
distribution from dAu interactions. 

6. As it is found, see Fig. El the pseudorapidity density at midrapidity shows increase in 
the difference between the values obtained from central nucleus-nucleus collisions and those 
from pp/pp interactions as the cm. energy of reactions increases above the highest SPS 
energy. As shown above, this difference is connected with a size of the Lorentz-contracted 
volume. At the cm. energies available today in nucleus-nucleus experiments, this difference 
is comparable with the difference in the measurements. It is of interest to know the density 
values for higher than the experimentally reached a/snn i n order to check whether the devia- 
tion between the pp/pp and nucleus- nucleus midrapidity densities becomes more pronounced 
as it is expected. Unfortunately, no measurements, except the midrapidity density for pp/pp 
data, are available for the variables in Eq. (J2J beyond a/snn = 200 GeV. The densities for 
pp interactions are measured up to y/s^p~ =1.8 TeV, and these high-energy pp data can be 
compared to at least some estimates for the midrapidity density in nucleus-nucleus collisions 
at v /snn = y/s pp /3. To this end, we extrapolated the values calculated from Eq. (J2J) utilizing 
the function analogous to that found |Hj to fit well the pp data. Both the prediction for 
nucleus-nucleus collisions and the fit for pp data are shown in Fig. |21 

The dependence of the midrapidity density on the cm. energy in nucleus-nucleus in- 
teractions, which we show by the dashed line in Fig. El agrees well with the calculations 
we made above based on Eq. (J2J) shown by solid line. This dependence indicates that the 
midrapidity density in heavy ions increases faster with the cm. energy than that in pp/pp 
collisions considered at the three times larger cm. energy y/s pp . As we predict within our 
consideration and show in Fig. |2]in an assumption that the same behaviour as at the SPS- 
RHIC energies extends up to the LHC energies, the density Pnn(0) is found to be ~ 7.7 
at the cm. energy ^/snn = 5.5 TeV of heavy-ion collisions at LHC. From the CDF fit [3] 
and assuming it covers LHC energies, one obtains the value of ~ 6.1 for the pseudorapidity 



7 



density at midrapidity expected for pp interactions at LHC at */s pp = 14 TeV. The value 
we obtain for Pnn(O) at LHC energy is consistent with that of ~ 6.1 given in the PHENIX 
extrapolation [HH ITT] within 1-2 particle error acceptable in the calculations we made. Our 
result is in a good agreement with the recent ATLAS Monte Carlo tuned values [62 j . Taking 
into account that the LHC ^/snn is close enough to v /s^ p "/3, the nearness of the values, 
predicted for the LHC by us and estimated independently and based on the experimental 
fit in the wide range of energies by PHENIX |TIH ITT] , demonstrates that our interpretation 
provides an experimentally grounded description and predictive ability. It is interesting to 
note that as the energy changes significantly the particle production process seems not to 
change from that above highest SPS nuclear collision energies, in contrast to what happens 
when the a/snn departs to the SPS energies from lower ^/snn as discussed. 

As soon as we describe nucleus-nucleus midrapidity density using pp/pp data in wide 
energy range spanning form a few GeV to hundreds GeV, now we can solve Eq. (J2J) for the 
mean multiplicity iV c h/(0.5iVp art ) to describe the nuclear data where no measurements exist 
and to predict the mean multiplicity energy dependence at higher energies. In this calcu- 
lations, we use the fits of p pp (0) [5] and jl] and approximation for pnn(O) as functions 
of the cm. energy shown in Figs. [T]and|21 From the resulted curve for A^ h /(0.5A 7 part ) as 
a function of a/snn plotted in Fig. ^ one finds that the heavy-ion mean multiplicity per 
participant pair at a/snn = 5.5 TeV is just about 10% above the A^f(^/s pp ) approximation 
at y/s pp = 14 TeV and about 3.3 times larger the data mean multiplicity from heavy-ion 
collisions at ^/snn = 200 GeV. Again, this number is comparable with the estimate made 
from the pseudorapidity density spectra by PHENIX [HI1III] an d points out to no evidence 
for change of behaviour in the energy dependence as the ^/snn increases by about two mag- 
nitudes from the top SPS energy. Nevertheless, one can see that the data obtained at the 
highest RHIC energy give a hint to some border-like behaviour of the mean multiplicity 
where the pp/pp data saturate the nuclear data, and another region of the rise is possible 
to be found (as at low energies). This makes heavy- ion experiments at ^snn > 200 GeV of 
particular interest. 

Given the importance of the multiplicity and produced particle spectra as control observ- 
ables of the particle production process, various theoretical descriptions have been proposed 
and confronted to the nucleus-nucleus collision data, in particular to the mean multiplic- 
ity and the midrapidity density measurements |3llI2llISlinSl , the bulk variables considered 
here. Most of the models are shown to reproduce reasonably at least some part of the mea- 
surements or the trend obtained in the data, while they are found to be less successful in 
describing all the experimental findings. Despite many questions are still open and there 
are significant challenges for theory, the message of inability to reproduce the measurements 
taken at highest RHIC energies by scenarios solely based on interactions between hadronic 
objects is clearly given. The system produced at highest density and temperature ever 
reached experimentally requires a partonic approach to be explored for its description. In 
this sense, our interpretation invokes partonic picture being the necessary ingredient as the 
description utilizes the constituent quark framework. 

7. In summary, we analyse the similarities of the bulk observables obtained in heavy-ion 
reactions, (anti)proton-proton and e + e~ interactions in the large cm. energy interval from 
a few GeV to hundreds GeV. The similarities of such bulk observables like the mean mul- 
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tiplicity and pseudorapidity density at midrapidity are compared to our expectations under 
assumption of the universality of a mechanism of the multiparticle production in different 
types of interactions. Within the description proposed, secondary particles are produced ac- 
cording to the amount of energy deposited by participants and depending on their type. The 
bulk observables in headon nuclear collisions are treated on the base of interactions of nucle- 
ons, considered as effectively structureless patterns wholly depositing their energy, similar to 
that in e + e~ annihilation. It is shown that the observables obtained from nucleon-nucleon 
collisions are similar to those from nuclear or e + e~ interactions when one considers the single 
constituent quark pair interaction picture for nucleon-nucleon collisions at the cm. energy 
three times larger the energy of nuclear or e + e~ collisions. The correction to take into account 
the tripling factor in the cm. energy and for proper treatment of participating patterns is 
found to explain the actual measurements. Two separated regions in the measurements with 
a boundary about highest SPS nuclear experiment energy are reproduced. The approach 
is extended to and the consequences are made for particle-nucleus collisions showing agree- 
ment with recent RHIC data. The predictions for higher LHC energies within the proposed 
interpretation are made, found to follow the experimentally obtained extrapolations. As 
shown, experiments at RHIC, given a remarkable opportunity to expand investigations to 
the highest so far achieved density of matter, allow to reveal signatures of new phenomena 
and provide indispensable knowledge on the dynamics of strong interactions. This makes 
of great interest further comprehensive analyses of the data obtained up to RHIC energies 
and to be obtained at the higher energies challenging to the better understanding of the 
discoveries made and testing the predictions. 
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Figure 1: The charged particle mean multiplicity per participant pair as a function of 
the cm. energy. The solid and combined symbols show the multiplicity value of the most 
central heavy-ion (AA) collisions at RHIC as measured by PHOBOS Collaboration (■) in 
and by NA49 Collaboration at CERN SPS |2E] (*) and by E895 Collaboration at 
AGS [29J (H) (see also [2J, and, in pp collisions, the measurements made at CERN by UA5 
Collaboration (A for non-single diffractive, ▼ for inelastic events) at ^s pp = 546 GeV [371] 



and 



200 and 900 GeV [2T| and, at lower cm. energies, in pp collisions obtained at 



CERN-ISR (•) [32J and from bubble chamber experiments [331011 (B) , the latter compiled 
and analysed in [33]. (The inelastic UA5 data at -/s pp = 200 GeV is extrapolated in [T3] from 
the limited rapidity range to the full one.) The open symbols show the e + e~ measurements: 
the high-energy LEP mean multiplicities averaged here from the recent data values (o) at 
LEP1.5 = 130 GeV in jSEllSZI and LEP2 JsZ = 200 GeV in [32EBI, and the lower- 
energy data as measured by DELPHI [HHj (□), TASSO [H] (a), AMY jUJ (o), JADE g2] 
(+), LENA HH| (*), and MARK1 [Sj (#) Collaborations. (See also HSligEl for data on 
e + e~ and pp/pp collisions). The solid line shows the calculations from Eq. (J2J) based on our 
approach and using the corresponding fits (see text). The dashed and dotted lines show the 
fit to the pp/pp data from [30] and the 3NLO perturbative QCD fit [IB] to e + e~ data by 
ALEPH [3Z]. The arrows show the expectations for the LHC. 
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Figure 2: Pseudorapidity density of charged particles per participant pair at midrapidity as 
a function of the cm. energy of collision. The open and combined symbols show the pseu- 
dorapidity density values per participant pair vs. cm. energy per nucleon, a/snn, measured 
in the most central heavy-ion collisions at RHIC by BRAHMS 50, (o), PHENIX [10, (a), 
PHOBOS IHIEII (□), and STAR [HI] (*) Collaborations, and the density values recalculated 
in PU from the measurements taken at CERN SPS by CERES /NA45 [52 (+), NA49 [53] 
(0), NA50 [5Jj (o) and WA98 [55J (*) Collaborations and at Fermilab AGS by E802 and 
E917 Collaborations j5H] (*)• The nuclear data at ^/s^ around 20 GeV and the RHIC data 
at v /snn = 130 GeV and 200 GeV are given spread horizontally for clarity. The PHENIX 
data at a/snn = 62.4 GeV is taken from [H]- The solid symbols show the pseudorapidity 
density values vs. cm. energy y/s pp /3 as measured in non-single diffractive pp collisions 
by UA1 57J (■) and UA5 [1[3U] (*■) Collaborations at CERN SPS, by UA5 at CERN ISR 
(s/ s pp = 53 GeV), by CDF Collaboration at Fermilab |3] (▼), and in inelastic pp collisions 
from the ISR [22] (★) and bubble chamber [H^ loT?] (•) experiments. The data from the bub- 
ble chamber experiments |34^ I59*] are given as recalculated in [lj. The solid line connects the 
predictions from Eq. The dashed line gives the fit to the calculations using the 2nd order 
log-polynomial fit function analogous to that used [3] in pp data. The fit function from [5] 
is shown by the dashed-dotted line. The dotted line shows the linear log approximation of 
UA5 to inelastic events [I]. The arrows show the expectations for the LHC. Note that e + e~ 
data at y/s^ = 14 GeV to 200 GeV (not shown) follows the heavy- ion data [Sj. 
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